ESTIMATES AND IDENTITIES FOR THE AVERAGE 
DISTORTION OF A LINEAR TRANSFORMATION 

IGOR RIVIN 

Abstract. Let A be a linear transformation A : R" ^ ]R". We give 
sharp estimates on 



/ 

Js«= 



log||AM||dM 

, We also show asymptotic results (for large n) and evaluate a class 
of integrals over the sphere, including the integral of the logarithm 
of absolute value of one coordinate projection. 



Introduction 

Let A : R" — > R" be a non-singular linear transformation. In this 
note we estimate the "average distortion" of A. More precisely, we 
estimate 



(1) J(A) = -f log\\Au\\du. 

Js"-i 

We remind the reader that j- denotes the mean of the integrand with 
respect to the measure denoted by du, which, in our case, is the 
standard rotationally invariant measure on the sphere. Another way 
of putting it is that the measure is normalized to be a probability 
measure. This average is of considerable interest in dynamics (see 
the author's paper [IJ and references cited therein). A somewhat 
related problem of estimating the average 



(2) f \\Au\\du 



f 
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is one of the principal problems addressed in the paper Q, and we 
happily carry over some of the techniques and observations to the 
current setting. 

The main estimate is summarized in the following 

Theorem 1. Let the singular values of the matrix Abe Oi, . . . ,a„. Further, 
let 



H: N ^ IR 



be defined as follows: 
H(n) 

Then 



-log2-lZt\l neven, 

(n-l)ll _J_ 
k=\ 2k-l 



^-Lfjl^'^jh nodd. 



log n > £^ ^ log \\Au\\du - i log ^ j > 

This will be shown by way of Theorem |2l in Section [l] followed by 
the explicit computation of the lower bound in Sections |2l and |3l 

In Section m we will indicate asymptotic results (a "law of large 
numbers") which indicates that the upper bound in Theorem [T] is a 
better guess for reasonably well-conditioned matrices A. Note, how- 
ever, that the difference between the lower and upper bounds is 
asymptotic to y -I- log 2, and so the gap between the two bounds is 

dimension independent, which indicates that 5 log ^X!,r=i '^f) is the 
"right" approximation to I{A). 

1. A SHARP INEQUALITY 

Our first observation is that we can assume that the matrix A (as 
in eq. ^) can be assumed to be diagonal, since the average in eq. ^ 
does not change if the matrix A is replaced by the diagonal matrix of 
its singular values Ci, . . . , On- (see 0111 for more discussion of this), 
so we can replace the integral to be estimated by 



h 

We will estimate the quantity 



du. 



n 

2 



i=l 



AVERAGE DISTORTION 



3 



Since /e is scale invariant, it will be enough to estimate it (or, what is 
the same, Iz) on the simplex H: Yj1=\ o]- Let us change variables, so 
that Si = cr?. By the concavity of the logarithm function, the integrand 
of Iz is concave (in the (si, . . . ,s„) variables), and hence so is Iz itself. 
Since, in addition, z is symmetric, it follows that the maximum of Iz on 
H is attained when Si = • • • = s„ = 1 /n, while the minimum is attained 
at (any) vertex of H, for example when Si = 1, whle Si = 0, ii^l. 
We have just proved the following 

Theorem 2. 

log n > ^ log I^J, o]x] du-^ log I^J, ^ ^ log \xi\du. 

2. How TO INTEGRATE OVER THE SPHERE 

Theorem 3. Let / : K." ^ K. have the property that 

(3) f{ax) = g{a)+f{x). 
Let [I be the standard measure on S""^ Then 

(4) 2"/2-ir(n/2) r f{x)d^ = 

Js"-i 

J" /(xi, ...,Xn) exp I — i — ^ j dxi... dx„-cOn-i g{ry-^e-'^'^dr, 

where 

is the area o/S""^. 

Proof. Let us evaluate the first integral on the right hand side. First, 
let us transform to polar coordinates: 

(5) I = f{xi, ...,Xn) exp I — - — ^ j dxi . . . dXfi = 

I r"-^e-'''^dr ( f{rx)d^ = 
Jo Js"-i 

I r"-V'/2dr I (/(:c) + ^(r))d;u = 
Jo Js«-i 

I ^n-ig-r2/2^^ I f{x) + cOn-i \ r"-^ g{r)e-'''^dr, 
Jo Js"-i Jo 
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where we have used Eq. ©, and have denoted the area of the unit 
sphere S""^ by (jJ„-i. To evaluate oj„-i (and, at the same time, the first 
integral in the right line) let /(x) = 1. Then g{r) = 0, and we have the 
equation: 



r (-x^-...-xl\ r°° 

I exp dxi . . . dx„ = co„-i ( r" 

Jr" \ 2 / Jo 

The left hand side equals 

.oo \n 



(£ 



The integral on the right hand side can be evaluated by changing 
variables to u = thus getting: 

Jo Jo 
Finally obtaining: 

CO„-i - 



□ 



3. An extended example 

Let /(xi, . . . = log(lxil)). Then f{a{x)) = log(fl) + /((x)), so g{r) = 
log(r). Let us evaluate the integrals on the right hand side of (HI). First, 

(6) r log(|xi|)e-(^i+-+^">/2dxi ...dx„ = 
Jr" 

g-(x?+...+.2_^)/2^^^ . . . dx„_i ( log(|x|)e-'^'/2rfx = 

log(x)e-^'/2dx = 



-2"/2-i7i"/2(y + log(2)), 

where y is Euler's constant. The second integral is also simply eval- 
uated in terms of special functions: 



f 

Jo 



lo§{r)r"-^e-''^ = 2"'^-^T{n/2){lo§{2) + ip{n/2)), 
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where \p{x) is the logarithmic derivative of the T function. Putting 
everything together, we get: 



(7) 



2«/2-ir(f) 

_2(n/2-i)^n/2(^ + log(2)) - ^2-/^-^1 (|) (log(2) + ^(n/2)) 



7i"/2 n"'^ 



(y + log(2)) - — -(log(2) + i/^(n/2)). 

r(f) r(f) 

It is sometimes more useful to compute the mean of a function over 
the sphere, and we can do that too: 



(8) 



Js"-i 



^ (-(y + log(2)) - (log(2) + ipin/2))) 



i(-21og2-y-i/^(n/2)). 

Lest the reader is discomfitted by the appearance of the digamma 
function \p, we note the following simple formula for its special values 
at integer and half-integer points: 



_l-y + j:"iT'l' "even. 



4>{n/2) = 



-y-2\og{2) + 2Zti^'^^^, otherwise 



Curiously, this indicates the the mean of the log |xi | over the unit 
sphere is either rational or in Q[log(2)], since y always cancels.^ 

4. Laws of large numbers 

The methods of Q go through without change to show the follow- 
ing results: 



Good thing, since it is not known whether Euler's constant is rational. 
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Theorem 4. Let o-i, . . . ,o„, . . . be a sequence of -positive numbers such that 

lim = 0. 

Let An e GL(n, R) be a matrix with singular values 0\,..., o„.Then 

lim £ ^ \\Au\\du - i log|X j + ^ logn = 0. 

Corollary 5. The conclusion ofTheorem^holds under the assumption that 
there exists a constant c such that Oi/Oj < cfor all pairs i ^ 

Remark 6. The hypothesis of Corollaryl^says that the condition num- 
bers of the matrices A„ is uniformly bounded. 
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